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+^ ! Abstract 

A generalization of Termo Field Dynamics to a curved geometry is proposed. 
\Q • In particular a neutral scalar field minimally coupled to gravity is considered as 
matter content in a Robertson- Walker metric. A non linear amplification in the 
particle creation is obtained, due to the altogether action of thermal and geometric 
' effects. As a consequence the frequencies in the system look like red shifted with 
respect to the case where the thermal creation is not taken into account. 
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■ 1. Introduction 

| A realistic description of the dynamics of our universe must take into account 

all the possible interactions. But that is practically impossible. An approach that 
we can used is to separate the universe in system and enviroment. Then we supouse 
degrees of freedom with different hierarchy, some of them considered as quantum 
degrees of freedom and the other ones constitute the classical background. The 
interaction between both kinds of degrees of freedom can be produced through the 
changes in the spacetime which are driven by the dynamical Einstein equations 
and by spontaneous thermal fluctuations of the enviroment which acts as a thermal 
reservoir. In the cosmological semiclassical approach (see ref. [1]) the role of the 
classical background is played by the gravitational field. This background can also 
play the role of thermal reservoir, as is shown in ref. [2] and the conformal quantum 
thermal fluctuation of the gravitational field can be interpreted as the tilde fields 
of the formalism of Thermo Field Dynamics (TFD) [3]. However the tilde fields 
"inhabit" also in a curved geometry, then they are also affected by the changes in 
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the metric. So we will consider a matter field in a thermal bath, due to the radiation 
background, and also to the interaction with the gravitational field. Therefore we 
can use Quantum Field Theory in Curved Space Time at Finite Temperature , or 
in a more realistic way, instead of finite temperature, a nonequilibrium description. 
In order to do that it is necessary to introduce some kind of interaction between 
the system and the thermal bath. A simple and elegant model of that interaction is 
proposed in the TFD formalism [3]. The fundamental hypothesis is: the increase of 
the energy in the system can be produced by the excitation of an additional quanta or 
by annihilation of holes of particles from the reservoir. This is spontaneous process, 
as is stressed in ref. [4]. The other mechanism that increase the energy of the 
system is due to the creation of particles from the interaction between the scalar and 
gravitational fields (see ref. [1]). This interaction is introduced indirectly by means 
of the curved spacetime (CST). The evolution of the particle creation due to the last 
process is driven by the field equation. When the in-out situation is considered the 
Bogoliubov transformation that connect both Cauchy surfaces produces squeezed 
states as is shown in references [5] and [6]. In that sense the isotropic universe, 
taking this work into account, is another example, as the ones given in ref. [2], that 
relate the dynamics to the squeezed states. 

The TFD hypothesis, as is shown in ref. [3] is implemented introducing an 
extended Fock space that includes the tilde states \n > of the reservoir with the 
states | to > of the system. Then we have the extended Fock space given by: 

{\n,h >} = {\n >} <g> {\n >} 

In particular we use for the vacuum the notation |0 >= |0,0 >= |0 > |0 >, 
with the operators ak, ak, a^ and aj., such that 

a k |0 >= 0, a k |0 >= 

Ok^|0 >= | lk >, etc ajjO >= |lk >, etc. 

Moreover the operators satisfy the commutation relations 

[a k , ak' 1 '] = 5k,k', [5k, &k'] = ^k,k' 

In the formulation given in ref. [7] the tilde field is a fictitious one, it does not 
represent any physical magnitude. In ref. [8] Israel found a parallelism between this 
approach and the problems where an event horizont is present, like for example the 
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Rindler observer and the black hole radiation. Then the role of tilde fields is played 
by physical but hidden modes. In ref. [2] the conformal fluctuationof the metric 
was proposed as tilde field . The real role of the tilde field surely must be played by 
a more complicated combination of the different fields present in the background. 
The use of one particular candidate can be interpreted as the choice a kind of coarse 
graining in the model of the system in a thermal bath. In the present work we can 
obtain interesting conclusions without specifying the nature of the tilde fields. In 
order to calculate statistical mean values of the physical observables, in TFD (see 
ref. [3]), the thermal vacuum state is introduced. If T = T(f3) is the temperature 
of the bath, the thermal state is defined by the transformation 

\0,P>=U[a k , a k , a k f , a|]l° > (1.1) 

The operator U allows us to obtain the thermal vacuum as a mixture of states 
at zero temperature. Then the fundamental hypothesis of TFD can be imposed by 
the expression: 

a\O,0 >= Aa f |0,/3 > (1.2) 

In general for a system out of equilibrium, the operators are time dependent 
and related with the latest ones by transformations as follows (see ref. [9]): 

a(t) = S'^^aSit) 

aP(t) = S- 1 (t)a!S(t) 

Where the ft symbol means that in general (a"y ^ a because S is not neces- 
sarily unitary. The eqs (1.2) turn therefore to 

a(t)\O,0(t) >= X(t)a^(t)\0,/3(t) > (1.3a) 
< 0, 0(t) |a tt (t) = 7 (£) < 0, P(t)\a(t) (1.36) 

We can use in the following the notation of ref. [9], i.e. : A = F, 7 = fF~ x . 

In TFD the statistical mean value of some physical observable is the thermal 
vacuum expectation value of the operator associated to the observable (see ref. [6]). 
In particular the mean value of created particles will be 

n tp (t)=<0,/3\a^(t)a(t)\0,l3> (1.4) 
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The eqs (1.3) motivate the definition of the operators a{(3) and a* {(3) such that 



a(p)\0,p >= 
< 0,P\a!(J3) = 



(1.5o) 
(1.56) 



with a(P(t)) oc (a(t) - Fa tt (t)) and a*(P(t)) oc (a tt (t) - fF' 1 ). Following refs 
[2] and [4] we can choose / and the proportionality constant such that: 



1 —F 

-fF- 1 1 



a(t) 

ait(t) 



(1.6) 



/, F and are also time dependent functions. In order to calculate ntp by 
means of eq.(1.4), we need perform the inverse operation to transformation (1.6), 
i.e.: 



^(t)\_ 1/2 ,( 1 F\fa(p(t))\ 

~aHt))- {1 + n ^ /{1 - f) [fF-i lJU(W 



(1.7) 



2.Quadrivectorial notation 

It will be convenient to introduce a quadrivectorial notation, in order to include 
also the and a operators. Then we can define (see ref. [10]): 



A k := 



/ a k \ 
a k f 
a k 
\al/ 



(2.1) 



Let us also introduce the 4x4 matrix Y so that the transformation can be 
written as 



A k (/?(£)) = Y k (/3(t))A k (t) 



(2.2a) 



A k (t) = Y k " 1 (/3(t))A k (/3(t)) (2.26) 

where 

Y k (/3(t)) = (l + M 1 /^I ^ (23) 
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r^ = (l + n t e)-i/\l-f)-i(* L " L ) (2.4) 

Moreover I is the 2x2 identity matrix. 

From the aplication of the transformation (1.7) in eq. (1.4), we obtain 

n (l + n )=f(l-f)- 2 



then the relation between n t p and / is: 

ntp = (2.5) 

It is interesting to note two things: 

a) nt/3 does not depend on F. 

b) The Planckian spectrum corresponds to a particular / function; / = exp(—(3e) 
with e the energy by mode. 

3. Dynamical and thermal effects in CST 

The change in time of the operators is driven by the dynamical equations. 
These equations, in our case, has the form of Klein-Gordon ones in a curved space 
time. The main difference with flat space time is that the coefficients in a normal 
mode expansion are time dependent. In order to obtain the time dependence of 
the annhilation-creation operators we will use the method developed by Parker [11] 
about which we will present a brief review. In our case the system and reservoir is 
considered in a curved geometry. The system is for us a massless and neutral scalar 
field minimally coupled to the gravitational field, given by a Robertson- Walker 
metric. 

The action for the matter field is 

S=\j ^9d A xd^<p (3.1) 

and the metric 

ds 2 = dt 2 -a(t) 2 (dx 2 + dy 2 + dz 2 ) (3.2) 
then the variation of the action gives the field equation 



V M c"V = 



(3.3) 



(with fi = 0,1,2,3). 

In order to make the calculation easier, we can use a discretization in the same 
form as in [11], i.e. we can introduce the periodic boundary condition ^(x+nL, t) = 
</?(x, t), where n is a vector with integer Cartesian components and L a length which 
goes to infinity at the end of the calculation. Then the integral J d 3 k can be replaced 
by (27r/L 3 ) ^ k . In order to expand a general solution of the field equation we can 
introduce the set of functions {(f)^(x)} U{4>* in( x )} defined by 

Mx)= M ) P>^ ^ x -fj {k - t ' ) '* ) - (3 ' 4) 

where W, a real function of k = |k| and t, has to be determined by the field 
equations and the boundary conditions. Then we can expand the field as 

p(x, t) = Y}Mt)Mx) + aV^kW] = ^M>k(aO + a^\(x)] (3.5) 

k k 

with ak := a^(t = t\) when t > t\ > to and with (as in ref. [11]) 

V'k(^) = h(k, t) expik.x (3.6) 



MM) = - /2 [a(M)*e J *° + /3(k,t)*e J *o ] (3.7) 
(La(t)) ' y2W 

a and (3 are kwown as Bogoliubov coefficients. In order the normalization 
condition, holds for the base of field equation solutions, those coefficients must 
satisfy 

\a\ 2 -\(3\ 2 = l (3.8) 
Eq. (3.8) is equivalent to the parametrization of the coefficients in the form: 

a(k, t) = e" l7 « (k ' t} cosh 0(k, t) (3.9a) 
(3(k,t) = e^ (k ' t) sinh6>(k,t) (3.96) 

Replacing eqs (3.5)-(3.9) in the field eq. (3.3) we have the system of eqs 

(1 + cosrtanh#)M + 2W7 Q = (3.10a) 
6 



MsmT + 2W6 = (3.106) 
where we define r = 7a + 7/3 - 2 / Wdt', and 



M = -\{W/W) + \{W/Wf - \{a/af - \{a/a) + ^~W 2 (3.11) 



From eqs (3.5)-(3.11) we can obtain the Bogoliubov transformation between 
the operators given by 

a k (£) = e *7«(k,t) cosh6>(k,t)a k + e ll ^ Kt) sinh0(k, t)a f _ k 



a t k (t) = e -^ (k ' t) sinh6'(k,t)a k + e- i7Q(k ' t) cosh6'(k,t)a t _ k (3.12) 

Then we can define the mean value of created particles between the Cauchy 
surfaces labeled by the times and t, in the form 

n t , (t)=<0\a^(t)a k (t)\0> (3.13) 

Using the inverse transformation of eqs (3.12) it is easier to see that 

n , t (t)=<0,t\a^a k \0,t> (3.14) 

Moreover from eqs (3.10) and (3.13) we obtain that n can be written in the 
form no,t = (g — 1) _1 with g a function of the particle model used, i.e.: 

cos 2 T 

The function g plays a role analogous to that of / _1 from eq. (2.5). For some 
particle model eq. (3.15) permits us to obtain Planckian spectrum. 

In the quadrivectorial notation the transformation given by eqs (3.12) is 

A k (t) = O k (t)A k (3.16) 



with 



n = 



( e^ k >*)cosh0(k,t), e^ k >*)sinh0(k,t)P 
" Ve-^( k '*Wh0(k,t)P, e-^( k ' t )cosh^(k,t) 

P is a parity operator, such that 

Pa k = a_ k 



and is the 2x2 zero matrix. The transformation given by the matrix O is 
unitary (see ref. [11]). Then in our case ft" = t because the dissipation is not intro- 
duced. Moreover as the dynamical change of the operators is supposed independent 
of the thermal one, so we have 

A(£, 0(t)) = Y((3(t))Q(t)A := A(t)A (3.17) 



therefore 



A=(l + n ,) 1 / 2 (l + tho) 1/2 (? I) 



with 



R 



e ii a e^tanhflP 
e-^tanh^P e~^ a 



-Fe-^tanh^P -Fe~ % ^ 
—fF~ 1 e ila - fF' 1 e^nanhOP 



T = 

By means of the last transformation we can reobtain nt t p in the form 

n t ,p=<0\a^(t,p)a(t,p)\0> (3.18) 

It is easy to prove that also the following equality is valid: 

< O t t,/3\aia\0,t,/3 >=< 0|a f (t,/3)a(t,/3)|0 > (3.19) 

Performing the calculation of n fjj g, using the transformation (3.17), results (as 
it is shown in ref. [10]): 

n t p = n t o + n 0f 3 + 2n to n /3 (3.20) 

As is mentioned in ref. [10], the increment in the number of created particles, due 
to the temperature is equal to the one produced in curved space at zero temperature, 
when there is a particle distribution in the initial state, as is shown in ref. [11]. 
That is due to the analogous role of the operators a k to the a_ k in curved space 
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(see ref.[12]). A similar amplification effect in the particle creation was obtained 
in ref.[13] with a massive scalar field interacting with a massless scalar field in a 
Robertson- Walker metric. 

4. Thermal equilibrium 

We will see now what the effect in the thermal distribution of the field modes 
is, due to the presence of the thermal reservoir. In order to do that the equilibrium 
situation is considered. Then let us extremize the Helmholtz free energy 

F=E-^K (4.1) 

Where E is the statistical mean energy at time t and temperature T = T{0), 
given by 

E=<O,t,0\H(t)\O,t,0> (4.2) 
whith H the hamiltonian operator. K is the entropy, which is obtained as 

K =< O,t,0\K(t)\O,t,0 > (4.3) 

where the K operator was introduced in ref.[7] in the form 

K = -^{a^a^log n tf3 (k) - a^a\log (I + n tp (k))} (4.4) 

k 

Therefore we obtain the entropy of a Bose gas (see ref.[I4]), given by 

K = ~ ^{ntpWog n tp (k) - (I + n tf3 (k))log (1 + n t0 (k))} (4.5) 

k 

with n t /3 as in eq.(1.4). 

We can calculate the energy, as in ref. [9] , from the tensorial energy- momentum 
operator 

= y^lakQk'-D^ [ipkjx) , ifa' {%)] + a k a t k ' J D M ^[V'k(a:),^*k'(^)] + 
kk' 

a\avD^[iP\{x),iPv{x)} + a\aK> D ^\{x) {x)\ + {k ^ k'}} (4.6) 

with _ ^ 

D^^p^] = -d^d^ip - -g^ipd^ 

9 



Then we can calculate the hamiltonian by 

H = I a s d 3 xT bo 



where f d 3 x = L 3 in the discrete approximation. 

Instead of eq.(4.2) we can calculate the energy, in an equivalent and more easier 
way, as E =< 0, /3\H[a k (t), a k t(t)]|0, (3 >, with 

H = ^K(t)a k ^(t)F 1 + a t k (t)a t k ,(t)F 1 * + a k (t)a t k/ (t)F 2 + a t k (t)a k/ (t)F 2 *} (4.7) 
kk' 

where 

1 • • Ik 2 
Fi = -0 k 0k' + — )0k0k' (4.8a) 



with = 0[W] given by eq. (3.4) and h the Hubble coefficient that for R-W 
universes is the relative velocity of the scale factor, i. e., h = a /a. If we suppose 
that W satisfies the Cauchy data that diagonalized the hamiltonian at t time (see 
ref.[15]), the expression (4.7) turns to 

H = Y,W{l + a k \t)av(t)} (4.9) 
kk' 

Taking into account eqs. (1.4) and (4.1) we have 
i r = ^W^ + ^) + ^[n^(k)^n^(k)-(l + n^(k))^(l + n t/3 (k))]} (4.10) 

In equilibrium 

SF = (4.11) 



Then the equilibrium condition is 

8F 



dntp 

That gives us the usual Planckian distribution 



(4.12) 



10 



If we call Wq the frequency of the normal modes in thermal equilibrium and 
curved spacetime, when the thermal effect, in the particle creation, is not considered, 
we can write in equilibrium n t Q = l/(e^ w ° — 1). Moreover from eq.(3.20) we have 

n tf 3 = n t0 (l + — + 2n 0( 3) 

For high temperature, when the approximation Wf3 « 1 (idem for Wq ) is 
valid, taking until the second term in the exponential expansions, gives 

Wn 

W ~ ^ (4.14) 

In the denominator of eq. (4.14) all the terms are positive, therefore the fre- 
quency W looks like red shifted with respect to Wq. 

5. Conclusions and Comments 

The particle creation due to the thermal bath and the one coming from the 
interaction with the geometry, both mechanisms of particle creation as we can see 
from eq. (3.20), act symmetricaly, however the nature of both processes is different. 
One of them is a spontaneous process due to the thermal fluctuations of the reser- 
voir, while the other one is produced by the dynamical change of geometry. The 
symmetry is due to the fact that in the two situations the vacuum state changes 
by means of Bogoliubov transformations, because still for the fields in the curved 
geometry the field equation is Klein-Gordon like. The similarity of the two trans- 
formations was analized in ref. [16]. There we have seen that the isotropy of the 
spacetime produces a mirror symmetry between the modes of the field and one of 
these symmetrical spaces can be associated with the tilde field of TFD. In a more 
complex case, when there is interaction with other fields, additional amplification of 
the particle production can be present (see ref. [13]), producing a "cascade" effect. 

From eq.(4.14) we can see that the observed frequency of the system with 
reservoir, with respect to the one at zero temperature looks like red shifted due to 
the effect of thermal creation of particles (or modes of the field). 

The calculation is identical if the field is massive, but the interpretation can be 
different, because the thermal effect changes the redistribution of modes of the field 
by the interaction with the reservoir. Then the particles obtain quanta of energy of 
the medium, but no new particle appears. In the massless case the simple creation 
of modes can be interpreted as new massless particles. In a flat spacetime it is 
necessary to introduce some symmetry breaking mechanism, in order to produce 
massive particles from the vacuum (see ref. [17]). However in curved spacetime 
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massive particles are created without an explicit breaking of symmetry. This can 
be interpreted using the analogy between the vacua in the curved spacetime and the 
ones related to accelerated observers in flat spacetime (it is known in the literature 
[1],[18], as Rindler observer). As it is shown in ref. [19] the Rindler observer has a 
Planckian spectrum in his reference system while in the Minkowski vacuum particles 
are not observed. It is similar to the behaviour of the electromagnetic field with 
respect to different reference systems. Actually we have treated with a toy model 
of electromagnetic the field, represented in our case by the massless scalar field. 
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